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Abstract 


We apply a new method to solve fractional partial differential equations 
(FPDEs) with proportional delays. The method is based on expanding the 
unknown solution of FPDEs with proportional delays by the basis of Bern- 
stein polynomials with unknown control points and uses operational matrices 
with the least-squares method to convert the FPDEs with proportional de- 
lays to an algebraic system in terms of Bernstein coefficients (control points) 
approximating the solution of FPDEs. We use the Caputo derivatives of de- 
gree 0 < a < 1 as the fractional derivatives in our work. The main advantage 
of using this technique is that the method can easily be employed to a variety 
of FPDEs with or without proportional delays, and also the method offers a 
very simple and flexible framework for direct approximating of the solution 
of FPDEs with proportional delays. The convergence analysis of the present 
method is discussed. We show the effectiveness and superiority of the method 
by comparing the results obtained by our method with the results of some 
available methods in two numerical examples. 
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1 Introduction 


Fractional partial differential equations (FPDEs) plays an important role in 
recent years. Many significant phenomena in mechanics, biology, networks, 
signal processing, reaction process, and other fields in sciences and engineer- 
ing are modeled by FPDEs. Although FPDEs are exploited to describe these 
phenomena carefully, but it is very challenging and irritating to obtain the 
exact solutions or at last approximate solutions of FPDEs. In recent years, 
many numerical methods are used to solve FPDEs. For sample, Cheby- 
shev collocation method [1, 6, 19], Legendre collocation method [29], Jacobi 
collocation method [9, 28], variational iteration method (VIM) [13, 26], ho- 
motopy perturbation method (HPM) [15, 23], Laplace transform method [4], 
Adomian decomposition method [25], and finite difference method [30]. 


FPDE with proportional delays is a particular type of FPDEs. This 
case of FPDEs models various phenomena such as shallow-water waves or 
nonlinear waves, which come from quantum field theory, medicine, biology, 
and climate models. For more details and various applications, we refer the 
interested reader to [20, 27]. Recently, some scholars have attempted to 
solve these types of equations. Sakar, Uludag, and Erdogan [20] applied the 
homotopy perturbation method to obtain an approximate solution of FPDEs 
with proportional delays. Singh and Kumar [22] used VIM for solving FPDEs 
with proportional delays. Shah et al. [21] suggested the natural transform 
decomposition method for constructing approximate solutions of FPDEs with 
proportional delays. 


Bezier functions are appropriate tools for obtaining the solution to these 
cases of problems. Easy definition, well performance, and interpolation of 
end-points are some virtues of these tools. Many scholars have applied con- 
trol points of the Bernstein—Bezier form in their papers. Zheng, Sederberg, 
and Johnson [31] introduced the Bernstein—Bezier form for solving ordinary 
differential equations. Chakrabarti and Martha [3] used Bezier curves for 
solving the Fredholm integral equations of the second kind. Ghomanjani and 
Khorram [5] used this approach for solving the quadratic Riccati differential 
equation. Karimi, Bahadorimehr, and Mansoorzadeh [8] proposed the use of 
this method to solve PDEs. Also, some researchers applied the Bezier control 
points method to obtain an approximate solution of fractional problems; one 
can see [2, 7, 10, 11, 14, 17, 18] and references therein. We suggest a tech- 
nique similar to the one used in [10], for solving two-dimensional fractional 
optimal control problems. 


In this paper, we consider the following FPDE with proportional delays 


Ou(a1@, byt) O™ Ulan, L, bn, t) 
Ox - Ox™ 


oDFu(x,t) = h(x, t, u(aox, bot), 


), @) 


with initial conditions 
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k k 
u) (x,0) = 8 (a), k=0,1,2,...,¢—-1, (2) 


where x,t € [0,1], ai,b; € [0,1] for 7 € {0,1,2,...,m1}, also s(x), k = 
0,1,2,...,q¢ — 1, where q € N, are smooth functions and h is the partial 
differential operator. Also the fractional derivative is defined in Caputo sense. 
The real numbers a,’s and 0;’s are given constant proportional delays. The 
necessary conditions for the existence of the solution of (1)—(2) and regularity 
properties of solutions can be found in [20, 21, 22]. 


Bezier control points with least-squares method are employed for solving 
FPDEs with proportional delays (1)—(2). The novelty of our strategy is in 
producing operational matrices —based on Bernstein polynomials— together 
with the least-squares method converting FPDEs with proportional delays to 
a simple optimization problem. 


This article is formed as follows: Section 2 introduces some properties 
of fractional derivatives and Bernstein polynomials. In Section 3, we define 
a class of FPDEs with proportional delays and employ operational matrices 
based on Bernstein polynomials to convert these equations to a mathematical 
programming problem in terms of some unknowns control points. In the 
fourth section, the convergence of the proposed method is studied. The 
applicability and efficiency of the proposed method are shown through two 
numerical examples, and comparison with some other methods is investigated 
in Section 5. Finally, the conclusion is presented in Section 6. 


2 Preliminaries 


This section contains some basic definitions and properties of fractional cal- 
culus, Bernstein polynomials, and function approximation. 


2.1 Some preliminaries in fractional derivatives 


Definition 1. Let w € R and let n € N. A real function f(x), 2 > 0 belongs 
to the space C,, if there exist a real number r > ys and g(x) € C[0, co] such 
that f(x) = x"g(x). Moreover, f € C7; if and only if fMEC,. 


Definition 2. The Riemann—Liouville fractional integral of order a > 0 
is defined as follows: 


of2(f(z)) = 7 [en 9ar 2 >0, 


Tia) 
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where f(x) € C, and I(a) = 5° r¢-Ye-Tdr is Euler’s gamma function. 


Definition 3. The Caputo fractional derivative of order a > 0 of the given 
function f(x) € C,, is defined as 


oDe f(t) = oly “D" f(x) 
1 


a : n-—a-1 f¢(n) = 
Teer (a — T) f(r) dr, n-l<a<n, 


where n is a nonnegative integer number; in fact n = [a]+1 while [a] indicates 
the integer part of a > 0. For a €N, the Caputo derivative operator is the 
usual derivative operator of integer order. 

We list some properties of fractional integral and fractional derivative for 
f€Cu, w> 1, and a, 8 > 0, as follows: 


oz (f(2)) = f(a), 
T(k+1) 


wl? (a*) = fetta” ke NUO0,2>0, 
oDF off (f(x) = f(), 
T(k+1 
oD? (2*) = ay kENUO,«¢>0,k > [a], 


n n 


0D (defile) = ae oDe fila), 


1=0 i=0 


oDz (C) =0, CER, 


where {c;}"_) are constants. 


2.2 Some preliminaries in Bernstein polynomials 


Definition 4. The Bernstein polynomials of degree n over the interval (0, 1] 
are defined as follows: 


Bre = (“Jao 
n! 


~ ino 
If one uses the binomial expansion for (1 — t)"~*, then 


ary => ("7') ("aye (3) 


k=0 


for i=0,1,2,...,7,, where (") 
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Definition 5. We define a Bernstein vector ¢,,(t) as 
dn(t) = [B3(t) BR(t) BR(t)... BRU]; 
the Bezier polynomial of degree n over the interval [0, 1] is defined as follows: 
Pa(t) = Cbn(t), (4) 
where 
C= [ts Gy Gye el, (5) 
is the vector of constant coefficients that we recall its entry as the control 


points. 
Thus by (4), 


Pb = \- GB): 
1=0 


Lemma 1. By using (3), we define 
bn(t) = dnTr(d), 


where T,(t)=[1 ¢ #2 ++» t™)? and vp = U"*? is an (n+ 1) x (n+ 1) 
upper triangular matrix that can be expressed by 


(—1)3-*n! ee 
n 7 a ee 
WT ee = 4 (n— J)! — i)! (6) 
0, t>J, 
where 7,7 = 0,1,...,n. 


Proof. See [2]. 


Theorem 1. Suppose that f(a) € C™[0,1]; then for each k € N, k < m, 
and k <n, there exists an (n+ 1) x (n+ 1) matrix D such that 


f(a) ~ P® (x) = CD*on(z), 


where C is shown in (5) and D = w,, AY, with 


Ace ade GH 941 
ee “i otherwise, 


for? =0,...,n and j = 0,...,n—1, and V can be expressed by 


Vy = 0553 k=1,2,...,n, 
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where 7, ;, is the kth row of #1. 


Proof. For details, see [8]. 


2.3 Function approximation 


Definition 6. For the product of two Bernstein polynomials of degree n,m € 
NU {0}, we define 
7 Bi" (a, t) = BP (ax) By (t), 
where t, x € [0, 1]. 
Lemma 2. The set 
Y = Span { 9 Bo" (#,t), 0 BY" (2, t),..-, 0Bm(#,t), TBo"(#,t),-.-, nBm (a, t)t 
is a complete basis for the Hilbert space L?[Q], where 2 = [0,1] x [0, 1]. 


Proof. Since Y is a finite subset of L?[Q], it is complete. For more details, 
see [12]. 


Lemma 3. If 
Y = Span{ § BG" (2, t), GBT (a, t),.--, 6B (at), TBO" (2, t),---, nBm (a, O)} 


and u(x,t) € L?[Q], then ¢2(x)Ud¢m(t) is the best approximation of u(z, t) 
out of Y, that is, 


u(x,t) Y Unm(2, t) Sa Has iat (a, t) a5) ue By (t) 


i=0 j=0 1=0 j=0 
= $7, (1)Udm(t), 
where U is the (n+ 1) x (m+ 1) matrix such that 
= Qn (¢n(2), (u(2,t), Pm(t))) Qn 
and Qn = UnGnvn', 


1 2 1 1 
2 3 1 
re ee! iia 
2 3 4 n+2 
Gn = r ‘ ’ 
1 1 ire 1 
n+1 n+2 n+3 2n4+1 


where 7%, is defined in (6). Also the inner product (-,-) on the space C(0, 1] 
is defined by 
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(f,9) = ‘i f(a)g(a)dx, fe), g(a) € C10, 


Proof. See [12]. 


To compute Caputo fractional derivatives for an arbitrary function f(z), 
we define an operational fractional derivative matrix of order a > 0 by using 
Bernstein polynomials. 


Theorem 2. Let ¢,(x) be a vector of Bernstein polynomials. Then for 
each a > 0, we have 


oDo on (x) = D bn(Z), 


while D° = wo, KOH, where K° is an (n+ 1) x (n+ 1) matrix, defined as 
follows: 


Keitl1jg+l= 


TG+1 
wo (B(pi, a) — BUSspisn)), 7ENUO, i> [a] or 
(i—a+1) t€NU0,7i>|lal,2=9, 
0, otherwise, 


where the floor function |a]| and the ceiling function [a] indicate the largest 
integer less than or equal to a and the smallest integer greater than or equal 
to a, respectively. Also, B(p,q) and B(x;p,q) are, respectively, the beta 
function and incomplete beta function, for i,j = 0,1,...,n, pj =i—[al], and 
qG = [a]—a+1. Also H® is the following (n+ 1) x (n+ 1) matrix: 


H° = |[H(1), H%(2),...,H%(n+1)], 
where H®(i+1) = (2'*1—-°, 6,(2)) A751, i =0,1,...,n, and An = (bn(2), bn(2)). 


The matrix D® is called the operational derivative matrix for the fractional 
Caputo derivative of order a. 


Proof. see [10] 


3 Approximate solution of FPDEs with proportional 
delays 


In this section, we consider the FPDE (1)-(2). One may assume that the ex- 
istence conditions for problem (1)—(2) are satisfied. We develop using Bern- 
stein polynomials for solving FPDE (1)-(2). To achieve this aim, by using 
the least-squares method, we convert problem (1)—(2) to a mathematical pro- 
gramming problem. For this purpose, define the residual function R(a,t) as 
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follows: 


Ou(a 2, byt) O™ Ulan, L, bn, t) 


a ae ). 
(7) 


R(x,t) = oDPu(a,t) — h(a,t, u(aoz, bot), 


It is clear that when we get the exact solution, then the residual function 
R(a,t) is zero for (a,t) € (0,1) x (0,1). To get the approximate solution of 
problem (1)—(2), we can follow Step 1 up to Step 6. 


Step 1: Choose natural numbers n and m and approximate the solution of 
u(x,t) by Bernstein polynomials. 
u(x,t) & Unm(x,t) = w(x, t)Gf (t)Udm(t) + &(2), (8) 


where w(x,t) and €(x) are given functions such that u(x, t) satisfies the 
initial conditions (2). If one chooses w(x, t) = t and €(x) = s(x), then 
we have 


Unm(2, t) = tol (x)Udbm(t) + E(2). 


But, by Lemma 3.5 in [2], the function t¢7(2)U¢m(t) can be rewritten 
as follows: 


tof (2)Udm(t) = $f (2) Udbm (tt = $F (e)Udbm (t) (Com(t))™ 


So, 


u(x,t) Y Unm(x,t) = 65 (r)Udbm(t) + €(z). (9) 


Step 2: According to the previous section, one can easily find that 
oDFu(az, t) ~ o§(a@)UD" dm (t), 
u(aox, bot) ~ 1 (apxz)Udm (bot) + €(aoz), 


Ou(a,2, bit) dt Ga DIG i+ acta): 
Ox dx 


O™ Ulan, X, bn, t) 
Ox 


a" (ani) 


= $3 (an,2)D"! Udm (bn t) + drm.” 


and 


ke 
u (0,0) © 6 (2D Udn(0) + EO = (a), k=0,1,....-1 
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Step 3: We express the residual function (7) by the Bernstein polynomials as 
follows: 


where the control points r;,; are functions of coefficients ujj;, 7% = 
0,1,2,...,n, 7 =0,1,2,...,m. In other word, 


Tig = Tig (uo,0, Uo,1; U0,2)+++) U0,m; U1,0;--- »Unym) : 


n m 
Step 4: Create a cost function F = >> >> r?,. It is obvious that the square 
i=0 j=0 
function F' is a function of uo,o, Uo,1, Uo,2,-- +; UO,m; U1,0;+-+;Unm- 


Step 5: Construct a mathematical programming problem as follows: 


minimize F = S- S- ri (Uo,0, U0,1; +++; U0,m,U1,0;++++Unym), (10) 
i=0 j=0 
s.t. 

Ss (Sram, } wrerBy 0 +e) = 40, 
i=0 j=0 \1=0 
k=0,1,2, ,q—1, 
O0<2<l, 
0<t<1 (11) 


Step 6: Optimization problem (10)—(11) can be solved by some useful computa- 
tional softwares such as MATLAB. By solving this optimization prob- 
lem, one can find the values of control points u;,; 1 =0,1,2,...,n,j7 = 
0,1,2,...,m. Then the approximate solution of u(#,t) could be found 
by using these control points in (8). We need to mention that we 
have solved optimization problem (10)—(11) by using discrete points 
Lp € (0,1) and ¢, € (0,1). 


4 Convergence analysis 


In this section, we show the convergence of the approximate solution that has 
been found by Bernstein polynomials for the FPDE (1)—(2) with proportional 
delays. In this way, without loss of generality, we suppose that n; = 1, that 
q = 1, and that 0 < a < 1 for the FPDE (1)-(2). Consider the following 
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problem: 


0 byt 
L |x, t, u(aox, bot), puoi it) 


) oDeu(z, t) 
Ou(a,x, bit) 


= oDfu(z, t) _ h(a, t, u(aox, bot), Ox ); 


(12) 
with initial conditions 
u(z,0) = s(x), (13) 


where (x,t) € 2 = [0,1] x [0,1], the constant proportional delays ao, a1, bo, 
and b; are in the interval [0,1], also s(x) is a given function and L is the 
partial differential operator. 


Consider (C'(Q), || - ||) as the normed space of all differentiable continuous 
functions, where || - || is defined as follows: 


Of 
Ox 


_an [tl 
I= Mo + |S] +] Zz» seem. 


Lemma 4. If = [0,1] x [0,1], then the defined function L : C(Q) —> R in 
(12) is uniformly continuous on (C(Q), || - ||). 


Proof. See [16]. 


Now, we prove the convergence of the approximate solution for the FPDE 
(1)-(2) when the degree of the Bernstein polynomials increase. For this in- 
tention, we remark the following lemma, that is an obvious result of the 
Stone—Weirestrass theorem. 


Lemma 5. If u(x,t) is a continuous real-valued function defined on Q and 
€ > 0, then there exists a sequence of Bernstein polynomials {unm (x, t)} 
and M €N such that 


nymeEN 


u(x, t) — Unm(2, t)|| < e, (z,t) EQ, nsm>M. 


Proof. See [24]. 


Theorem 3. Consider the functional L on C(Q). If qm is the minimum of 
the functional ZL on C(Q)(]Y, then 


TMnm <2? 0 


where LE and Y are defined in (12) and Lemma 2, respectively. 


Proof. Let u*(x,t) € C(Q) and € > 0 be such that L[u*(a,t)] < ¢. By the 
minimum properties, it is possible to find such a function u* (x,t). Consider- 
ing Lemma 4, if u(x,t) € C(Q) and ||u — u*|| < 6, then 
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|L{u] — L[u*]| <e. 


We use Lemma 5, for suitable values of n and m, there exist Unm(x,t) € 
C(Q) MY such that |Junm — u*|| < 6. If we show LlUnm] by nm, then we 
have 


Tnm = |L[Unm] — Llu*] + Llu"]| < |L[unm] — Llu"]| + |L[u"]] < 2e. 


Because € is arbitrary, this completes the proof. Note that the minimum 
value of the functional L in (12) on C(Q) is zero. 


5 Numerical examples 


In this section, we give two numerical examples and apply the method pre- 
sented in Section 3 for solving them. These test problems are solved by using 
powerful MATLAB 2018b software on an Intel Core i7-7500U. These test 
problems demonstrate the validity and efficiency of this technique. 


Example 1. Consider the following fractional partial generalized Burgers 
equation with proportional delay (see [20]): 


0? u(x, 5) t 


oDi u(z, t) = Ox2 ula, 5) u(x,t), 


where x,t € [0,1] ,0 <a <1, and the initial condition is 


u(x, 0) = 2. 


The exact solution of this example is u(x,t) = 27e' when a = 1. If one 
assume that t(2,t) is the approximate solution, then the error between exact 
and approximate solution can be shown by 


| u(x,t) — u(a, £) |, (a,t) € (0,1) x (0,1). 


The comparison of the absolute errors between exact and approximate solu- 
tion by employing HPM [20], AVIM [22], and the proposed method for n = 2 
and m = 5 is shown in Table 1. Figure 1 shows the graph of the approximate 
solution of u(«,t) for different values of a when x = 1. 
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Table 1: The absolute value of error for a = 1 in Example 1. 


x t | presented method HPM(20] AVIM([22] 
0.25 | 0.25 8.7365E — 08 5.30£ —07 | 2.1974E — 08 
0.25 | 0.5 5.6597E — 08 1.7735E — 05 | 1.4596E — 06 
0.25 | 0.75 3.6812E — 08 1.4087E — 04 | 1.7274 — 05 
0.25 | 1 1.7890E — 07 6.2178E — 04 | 1.0095E — 04 
0.5 | 0.25 1.0598E — 07 2.1230E — 06 | 8.7896E — 08 
0.5 | 0.5 1.4489E — 07 7.0943E — 05 | 5.8385E — 06 
0.5 | 0.75 6.2958E — 07 5.6348E — 04 | 6.9096E — 05 
0.5 1 1.3441E — 06 2.4871 E — 03 | 4.0379E — 04 
0.75 | 0.25 2.4535E — 07 4.7760E — 06 | 1.9777E — 07 
0.75 | 0.5 9.9618E — 07 1.5962E — 04 | 1.3137E — 05 
0.75 | 0.75 2.2682E — 06 1.2678E — 03 | 1.5547 — 04 
0.75 | 1 4.1328E — 06 5.5960E — 03 | 9.0853E — 04 


Figure 1: Approximated graphs of u(a,t) in a = 0.7,0.8,0.9,1 and 


Example 2. In this example, the following fractional generalized Burgers 


equation with proportional delay is considered (see [20]): 


such that 


oDfu(a, t) — 


?u(x,t) 


Ox? 


= 1 for Example 1. 
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where x,t € [0,1] andO<a< 1. 


The exact solution of this example is u(x,t) = ze’ when a = 1. The 
comparison of the absolute errors between exact and approximate solution 
by employing NTDM [21], HPM [20], AVIM [22], and the proposed method 
for n = 2 and m = 6 is shown in Table 2. Figure 3 shows the graph of the 
approximate solution of u(«,t) for different values of ~ when x = 1. 


Table 2: The absolute value of error for a = 1 in Example 2. 


x presented method | NTDM [21] HPM [20] AVIM[22] 
0.25 | 0.25 1.4976E — 06 2.1224F — 06 | 2.1230E — 06 | 8.7896 — 08 
0.25 | 0.5 7.1515E — 06 7.0943. — 05 | 7.0943E — 05 | 5.8386 — 06 
0.25 | 0.75 1.7332E — 05 5.6348.E — 04 | 5.6348F — 04 | 6.9096L — 05 
0.25 3.4298E — 05 2.4871E — 03 | 2.4871E — 03 | 4.0379 — 04 
0.5 | 0.25 1.0884E — 06 4.2550E — 06 | 4.24504 — 06 | 1.7579. — 07 
0.5 | 0.5 6.5797E — 06 1.4189 — 04 | 1.4189EF — 04 | 1.1677H — 05 
0.5 | 0.75 1.6607E — 05 1.1270F — 03 | 1.1270E — 03 | 1.3819% — 04 
0.5 3.3493E — 05 4.9742E — 03 | 4.9743E — 03 | 8.07584 — 04 
0.75 | 0.25 1.6021 — 06 6.3675E — 06 | 6.3670 — 06 | 2.6369. — 07 
0.75 | 0.5 7.5613E — 06 2.1283.E — 04 | 2.1283E — 04 | 1.7516£ — 05 
0.75 | 0.75 1.8304 — 05 1.6904E — 03 | 1.6905E — 03 | 2.0729. — 04 
0.75 3.6085E — 05 7.4614F — 03 | 7.4614F — 03 | 1.2114h — 03 


Figure 3: 


4.5 


Approximated graphs of u(x,t) in a = 0.7,0.8,0.9,1 and x 


= 1 for Example 2. 
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6 Conclusion 


The aim of this paper was to solve fractional partial differential equations with 
multiple proportional delays by using Bernstein polynomials. To reach this 
goal, we generated some operational matrices based on Bernstein polynomials 
for fractional-order derivative in the Caputo sense. Operational matrices, 
together with the least-squares method, were used in an algorithm for solving 
these kinds of problems. The steps of the algorithm were easy to implement. 
The convergence of the method was investigated. Two numerical examples 
were given to illustrate that this proposed method is accurate and precise. 
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